Buckled nano rod - a two state system and its dynamics 
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, We consider a suspended elastic rod under longitudinal compression. The compression can be used 

' to adjust potential energy for transverse displacements from harmonic to double well regime. The 

two minima in potential energy curve describe two possible buckled states at a particular strain. 
^ , Using transition state theory (TST) we have calculated the rate of conversion from one state to 

Q ■ other. If the strain e is between Sc and Asc, the saddle point is the straight rod. But for Sc < 4:6 c, 

the saddle is S-shaped. At Ec = iSc the simple TST rate diverges. We suggest methods to correct 
this divergence, both for classical and quantum calculations. We also find that zero point energy 
contributions can be quite large (as large as 10^) so that single mode calculations can lead to large 
errors in the rate. 
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^ , Considerable attention has recently been paid to two-state nano-mechanical systems Q, |M |^ Q H, IE Q and the 
(~| possibility of observing quantum effects in them. Roukes et al. P| proposed to use an electrostatically flexed cantilever 
, to explore the possibility of macroscopic quantum tunnelling in a nano-mechanical system. Carr et al. Q suggested 
• ■ using the two buckled states of a nanorod and investigated the possibility of observing quantum effects. A suspended 
^ elastic rod of rectangular cross section under longitudinal compression is considered. As the compressional strain is 
^ increased to the buckling instability 0, the frequency of the fundamental vibrational mode drops continuously to 
I zero. Beyond the instability, the system has a double well potential for the transverse motion (see Fig. ^1. The 
' ^ two minima in the potential energy curve describe the two possible buckled states at that particular strain Q and 
g , the system can change from one to the other, under thermal fluctuations or quantum tunneling. We use L, w and 
Q ■ d (satisfying L >> w >> d) to denote the length, width and thickness of the rod 0, 0, [^. F is the linear 
'— ' modulus (energy per unit length) of the rod and is related to the elastic modulus Q of the material by = Qwd. 

^ The bending moment k is given by — d^/12. We take the length of the uncompressed rod to be Lq. We apply 
^ , compression on the two ends, reducing the separation between the two to L. If y{x) denotes the displacement 
of the rod in the 'd' direction, then the length of the rod Ltotai — Jq dx^/l + {y'Y ~ L -\- 1/2 dx{y')^ . The 
compression causes a contribution to the potential energy VeiasUc = F/{2Lo) (Ltotai — Lq)^. In addition, bending of 
the rod in the 'd' direction cause bending energy Vb — Fk^ /2 dx{y"Y . Thus the total potential energy is given by 
^ : V[y{x)] = 1/2 dx{FK^{y"f + Fs{y'f) + F/{8L^){jl; dx{y'ff + F/{2L^){L - L^f . Here s = [L - Lo)/Lo is the 
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FIG. 1: Potential energy V as a function of the fundamental mode displacement Y. The shape of the potential energy is 
harmonic for e > Sc, quartic for e — Sc = critical strain [sc < 0) and a double well for e < Ec- 
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FIG. 2: The rod under compression: The central figure (A) shows the uncompressed rod of length Lq. On compressing to 
length L, the rod buckle, either to B- or to B+. 



strain. Extremisation of potential energy functional with respect to y{x) leads to 
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and the hinged end points have boundary conditions y{0) — y{L) = and y"(0) = y"{L) = 0. If e > Ec = — k^tt^/L^ 
then, the only solution to Eq. (Q) is y{x) — if e > Sc- But if e < Sc, two bucked states are possible. They are 
y{x) = ±Ay^2/L sin(7ra: /L), with A = ± y^2LoZ2(£7^^^)77r^. For e < Sc, all the normal modes of vibration about 
these are stable. The solution y{x) = is now a saddle point. One can calculate the barrier height for the process 



of going from one buckled state to the other over the saddle (linear geometry) as AEg 
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The kinetic energy of the rod is fi/2 ^y^dx, where /i = m/L^ is the mass per unit length. Using the boundary 
conditions for the hinged end points, we find the normal modes of the rod, y{x,t) = y„(a;)e*'^"*. At the saddle 
point, we obtain yn{x) = An \/'2/L sin(n7rx/L) , with n = 1,2,3.... The normal mode frequencies at the saddle point 



are given by ujI^^^ 



ujQ n — ^c, where ujq = tP'kI L'^ ^jFjJji. n = 1 is the unstable mode and it has the 



imaginary frequency i^^.^near i = i^Linear^ where VLnnear — \/ ^l^c — 1- For the buckled state, the normal modes 
are the same as at the saddle point, but the normal mode frequencies are different. They are = uJon^/v? — 1 
for n > 1, while wi = wo\/2 (e/ec — 1). The rate expression using classical TST is (Eq. 3.14 of reference [l3|) is 
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Ir/fc^ where N denotes the total number of transverse modes of 



the rod. One makes a negligible error by taking the value of N to be infinity and this leads to 
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The rate expression using quantum TST is given by 

( / T ■ '>l.~T \^ } 

j^quantum ^^Lmear^ 
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The above rates have a problem. As ^Je/ec 2, L02 ^ i) the second buckling instability sets in and the rates diverge. 
As the rod is compressed, first the mode Aiyj2/ Lsvuiirx / L) becomes unstable and this is the first buckling instability 
and the rod buckles as a result of this. The length at which this occurs shall be denoted by L/. If one supposes that the 
rod is compressed further keeping the straight rod configuration, then at a length L^, the mode A2\/2/L sin(27rx / L) 
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FIG. 3: Plot of AEg^^^lZr against (L - Lo) for a silicon rod of dimensions Lo = 1000 A, w = 200 A and d = 100 A. For this 
rod the first three instabilities occur at L/ — Lo = —8.364 A, Ls — Lo = —35.354 A and Lt — Lq — —89.238 A respectively. The 
first buckling instability is shown by an arrow, the second and third buckling instabilities are not shown in figure. 
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FIG. 4: Plot of AEgl^li^^ against (L - Lo), for a silicon rod of dimensions Lo = 500 A,w = 20 A, d = 10 A. For this rod the 
first three instabilities occur at Lf — Lq = —0.1646 A, Lg — Lo = —0.6597 A and Lt — Lq = —1.4893 A respectively and all 
three are shown in figure using crosses. Solid line is for linear saddle point (valid in the regime L > Ls), dashed line is for bent 
saddle point (valid in the regime L < La). 



too would become unstable and this is the second buckling instability. For e > Asc, there is only one saddle point. 

But for e < Asc, due to the second instability, the saddle point bifurcates into two. In a similar fashion one can have 
the third instability at a length Lt etc. In order to analyze the rate near and beyond the second buckling instability, 
we assume that the displacement has the form yo{x) = Aiy/2/Lsm{'Kx/L) + A2\/2/Lsm{2'Kx/L). Using this, the 
elastic potential energy is given by 

FnHAl+4Al)^ Fn^Al{e-e,) 2Fn^Al{s - 4s,) 
^^^1'^^) = mJLo + 2l5 + 15 • 

Finding the extrema of this potential leads to the following three solutions for {Ai, A2): (a) (0, 0): this is the 
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straight rod configuration. Between first and second buckling (i.e. Ls < L < Lf ), this is the saddle point. But 
after the second buckling, it is no longer a saddle, but it becomes a hill top. It has the energy Ehuitop = 0. (b) 
(±2/7r-\/LLo(ec ^ £)j 0): These are the buckled states and both of them have the same energy E}j = —FLo{£ — 
£c)^/2. (c) (0, ±l/7r-\/LLo(4ec — £))'■ These are the two new saddle points that arise from the bifurcation of the 
one that existed for Aec < s. At these saddle points, the rod has a bent (S shaped) geometry. These two have the 
same energy Ef^J^^i^ = —FLo(e — 4ec)^/2. Beyond the second buckling instability, the barrier height is given by 
^^Barrier — —3FLoec{—2e + 5ec)/2. Near the saddle, the normal mode frequencies are given by: (^g^nt i — *^Bent, 
^Bent 2 ~ ^o^/^is/^c ~ 4) and for n > 2, „ — LOQn\/ v? — 4. In the above ijj|jg„f i has an imaginary frequency 

with f^Berit — \/3i^o- Now the classical rate beyond the second buckling instability can be calculated taking the saddle 
to be the bent configuration: 



^classical ^ /££ £_^in^l^_2e+Be^) ^ (5) 
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It is interesting that the normal modes for this saddle retain their stability, irrespective of what the compression is. 
The quantum rate expression is given by 



Near the second buckling instability {^e/sc — > 2) vanishes, causing the rates in Eq. (|5J, Eq. Eq. |SJ and Eq. © 
to diverge. The cure for the divergence is simple for the classical rate and is given below in the equations (0) and ((HJ. 
All that one has to do is to include the quartic term in A2 of Eq. Q in the evaluation of partition function for the 
second mode at the saddle. All the other modes (at the saddle as well as at the reactant) are treated as harmonic. In 
the regime where L > Lg the rate is then given by (saddle is the straight rod) 




j^class^cal ^ (4 - ^ ) / ^^-li:! ( _ _ l)r(3 - . / -)r(3 + ./ -), (7) 

2L £c Z-T Lq£^ y TT/i £c \I £c V ^ 

where fint{o,, b) = dye . In the regime where L < Lg the rate is given by (saddle is the bent rod) 
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i?f™' = -r - l)\/z/»*(4- ^TTTTTl)- (8) 



Now we follow the work of Voth, Chandler and Miller for rate calculation using quantum mechanical transition 

state theory near the second buckling instability (where oj^inear2 small). The QTST rate under the harmonic 

approximation diverges at temperature T = hfl/{2kTr), but using this method it is possible to avoid the divergence 

UJ. For this we go beyond the harmonic approximations for the first two modes. The Hamiltonian for the first 

two modes (^1,^2) at the saddle may be written as: H = p\/[2y) +p2/{2^) + V{Ai,A2), where pi = fiAi is the 

.2 .2 

momentum operator canonically conjugate to At. The Lagrangian is L{Ai,A2) — ^/2[Ai + A2 ) — V{Ai, ^2)- The 
"centroid" partition function for these two coupled modes is given by the path integral Q* = J DA2{t) J DAi{t)6{Ai — 
Al)e~s[AAT),A2(rm, Here A denotes the position of the centroid, defined by A^ = 1/ {(3h) drA^^ir). Note 
that the centroid for the first mode is constrained at A\, where A\ is the value at the saddle [l3| and is equal to 
zero. We calculate the frequency of the unstable mode fl and the partition function for the second mode Q\ using 

a variational principle based on the trial action Striai = Jq^ dr {1/2 11 Ai (r) — l/2/i51^A^(T) + ^2(^2) + 

.2 

l/2/iyl2 (r) + 1/2iiuj2{A2){A2{t) — A2)^}. One can determine fi^, w|(^2) and ^2(^2) variationally, so as to get 

the best possible value for Q* . Once these values are obtained, we can proceed to calculate the rate, because in 

the trial action, the two modes are decoupled. Now the partition function of the second vibrational mode at the 

saddle Q2 may be approximated by Q2 = \J {fikT)/2Trh'^ J^^dA2e fe^, with an effective potential W2{A2) = 

A:Tlog[2fcT/(?ia;2(^2)) sinh(7iu;2(^2)/(2fcr))] + ^2(^2)- The tunneling current for the first mode may be taken as 
il/(27r sin(^^)), where Q is the variationally determined frequency of the reactive mode. In the regime where L > Lg 
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Classical Rate at T=300K 
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FIG. 5: Plot of rate of crossing from one buckled state to the other using classical transition state theory, for a silicon rod of 
dimensions Lq = lOOOA, w = 200 A, d = lOOA at T = SOOiS'. For this rod the first three instabilities occur at Lf-Lo = -8.364A, 
Ls — Lq = —35.354 A and Lt — Lo = —89.238 A respectively. The first buckling instability is shown by an arrow, the second 
and third buckling instabilities are not shown in figure. 



the rate may be calculated using (transition state is assumed to be straight rod) 

Rl-'- = "^"^tr:" sinh(|;i) sinh(||)Q; {[ "^^^^ ■ 
^sm(ig.) 2kT 2kT L% shM Tkr" ) 



In the regime where L < Lg the rate is given by (transition state is assumed to be bent rod) 

!e 2fcr huji , . , , huj2 
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= ^...^^) ^^"^(2fcr)""^(2fcr)Q^n (10) 
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Now we consider the Sihcon rods of different dimensions, as summarized in Table. ^ Si has an Young's modulus 
Q = ISOGPa and density p ~ 2230kg.m^^. First we consider a rod of dimensions lOQnrn x 20nm x lOnm, considered 
by Carr et al 0. In Fig. |31we plot the activation energy against the compression. At L ^ Lf, the activation energy 
is zero and as one compresses the rod, it increases rapidly, as it is proportional to (e — Ec)^- Over a very short range 
of compression (^ 0.14 A) it increases to about 6 kcal/mol. The classical rate is plotted in Fig. [3 At L = Lf, the 
potential is very flat and near this L the pre-factor (i.e., the attempt frequency) in the rate expression vanishes. This 
is the reason why the rate goes to zero as L ^ Lf. It is found that the rate increases at first as one compresses. 
This is due to the increase in the pre-factor for the rate from zero. Then the rate decreases due to the increase in 
the barrier height. For the quantum calculation, we choose the value of N to be equal to the transverse degrees of 
freedom that would be there if one considered an atomistic model for the rod. Thus, for this rod, we took N — 425. 
We also report (see Table Pi the quantum enhancement factor, defined by r=Quantum Rate/Classical Rate. The 
values of F was found to be « 1 implying that quantum effects are not important for this rod, within the observable 
range of compression. For a rod of dimensions 2000 nm x 20 nm x 10 nm, plots of the rates at 0.01 K is given in Fig. 
IHland it shows that quantum effects lead to an increase in the rate by about a factor of 10. However, the overall rate 
rapidly decreases to very low values as one compresses the rod by about 0.01 A (see the Fig. IHl and hence it would 
be very difficult to observe this quantum enhancement experimentally. We have also performed calculations for a 
hypothetical rod of dimensions 50 nm x 2 nm x 1 nm. It is yet not possible yet to have Si rod of these dimensions. 
However it should be possible to synthesize molecular rods of these dimensions m. In Fig. H we have plotted barrier 
height as a function of compression. In the regime L > Lg, the transition state has straight rod configuration and in 
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Classical Rate at T=300K 
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FIG. 6; Plot of rate of crossing from one buckled state to the other using classical TST, for a silicon rod of dimensions 
io = 500 A,w = 20 1, d = 10 A at T = 300 K. For this rod the first three instabilities occur at L/ - Lo = -0.1646 A, 
Ls — Lq — —0.6597 A and Lt — Lq — —1.4893 A respectively. Dotted line is for linear saddle point, dashed line is for bent saddle 
point and solid line include quartic terms for the second mode at saddle point. 



Comparison between Classical and Quantum Rate at T=0.01K 
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FIG. 7: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum TST, 
for a silicon rod of dimensions Lq — 500 A,w = 20 A, d = 10 A at T = 0.01 K. The solid line is the result using classical TST, 
dashed line is the result using quantum TST. 

the regime L < Lg, the transition state has bent configuration. Fig. shows the classical rate against compression, 
made at a temperature of 300 K. The rate obtained using the linear transition state is seen to diverge at the second 
buckling instability, but is finite for all L > Lg. Similarly, beyond the second buckling instability the rate is calculated 
using the bent saddle. Close to the instability, this rate too diverges, but a well behaved rate can be calculated using 
the approach of Voth et al. 0| outlined above. In the quantum rate calculation we have taken contributions from 
213 normal modes of this rod. We have compared the quantum rate with the classical rate at 0.01 K (Fig. d). There 
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Comparison between Classical and Quantum Rate at T=0.01K 
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FIG. 8: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum TST, 
for a silicon rod of dimensions Lo — 20000 A, w = 200 A and d = 100 A at T = 0.01 AT. The solid line is the result using classical 
TST, dashed line is the result using quantum TST. 



is a quantum enhancement in the rate of roughly 10^. This occurs as one compresses the rod by about 0.015 A. 
Again, fabricating such a rod and doing an experiment under such conditions is a formidable challenge. We now 
ask, what is the origin of a high quantum enhancement factor in some of the cases? This is not due to tunneling of 
the reactive mode [7 1 [al . but is a zero point energy effect. The earlier analysis took only the reactive mode in their 
calculations [slI^llOlllH. but our analysis takes all the modes. This leads surprisingly to a decrease in the effective 
activation energy. The effective quantum mechanical barrier height may be written as Aii'^™^"^^ + AEzero, where 
^Ezero is the difference in zero point energy of the modes between the buckled state and the transition state. For the 
linear transition state AEzero = J2n=2 ^i/n^ — e/sc — — 1^ = J2n=2 ^o{N — 1) (e/sc — 1) /2. Notice that 
AEzero is negative and hence leads to a lowering of the barrier height, which is proportional to N. This is the reason 
why the quantum effect is more pronounced for the 2000 nm bar (provided w and d are the same) as may be seen in 
Tabled The only possibility for observing quantum effects seems to be exciting the rod to a higher level in the buckled 
potential, suggested by Blencowe 7]. Thus, for the 2000 nm bar if one keeps L — Lq = —0.05 nm, the barrier height 
has the value 0.0935 kcal/mol. The classical rate at T = 0.01 K is 1.55 x 10~2038 ^g^-i j^^^^ Quantum enhancemenr 
factor r = 1.32 x 10^. Even with this enhancement, the net rate is far too small to be observed. But following the 
idea of Blencowe , the system which is initially in thermal equilibrium near the bottom of the well, can be excited 
to an energy level below the barrier maximum, then it will be more likely to tunnel through the barrier than being 
thermally activated over it. In such a case the rate analysis using only the reactive mode will be underestimating the 
rate by a factor of roughly 10^! 
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Dimensions 


Temperature 


Range of compression over which classical 
rate varies between 10^ and 10^ sec~^ 
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0.14 A 
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d=100 A 


O.OIK 
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300 K 
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lA 
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d=10 A 


0.01 X 


0.015 A 


« 10® 



TABLE I: Summary of results. 
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